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The problem of the lattice diffusion of two particles coupled by a contact repulsive interaction 
is solved by finding analytical expressions of the two-body probability characteristic function. The 
interaction induces anomalous drift with a vanishing velocity, the average coordinate of each particle 
growing at large times as t ly/2 . The leading term of the mean square dispersions displays normal 
diffusion, with a diffusion constant made smaller by the interaction by the non-trivial factor 1 — 1/n. 
Space continuous limit taken from the lattice calculations allows to establish connection with the 
standard problem of diffusion of a single fictitious particle constrained by a totally reflecting wall. 
Comparison between lattice and continuous results display marked differences for transient regimes, 
relevant with regards to high time resolution experiments, and in addition show that, due to slowly 
decreasing subdominant terms, lattice effects persist even at very large times. 

PACS numbers: 05.40+j, 66.30.-h, 71.35.-y 



I. INTRODUCTION 



Whereas classical ordinary diffusion of a single particle is universally known, diffusion of interacting particles 
does not seem to have drawn much attention. A notable exception is the seminal paper by Fisher JlJ, which introduces 
basic ideas and solves, in various cases, the problem of finding the probability for the reunion of a given number of 
"drunken walkers" wandering on a one-dimensional lattice. The so-called tracer problem also received some attention, 
following the solution given by Harris j|] for the Id case (recent bibliography on this subject can be found in Mallick' 
thesis H). More recently, intensive work Q, j^] and references therein, has been done one the asymmetric simple 
exclusion process ( ||) which describes the biased motion of a lattice gas with hard-core interaction. 

The interaction between diffusing particles appears to be relevant in many fields: one-dimensional hopping 
conductivity 0, ion transport in biological membranes ||, channelling in zeolithes ]l0|]. Generally speaking, 
the interactions are expected to play a dominant role in low-dimensionality systems and/or sytems with geometrical 
constraints. Such problems have been analyzed in the continuous space limit of the so-called sing le-ffle model 
& 

I here consider one of the simplest problems, namely that of two diffusing particles on a Id-lattice with a 
repulsive contact interaction, by directly solving the lattice master equation using elementary methods and obtaining 
the exact two-body probability at all times; this allows to find the behaviour Vi of the mean square dispersion of the 
coordinates. The totally asymmetric version of this problem (directed random walk in which each particle can move 
only in one direction) was treated in ref. H using a Bethe ansatz; among other results, these authors gave explicit 
asymptotic behaviour of the two first moments of the coordinates in the case of two diffusing particles. 

The continuous space limit of the model here considered is simply related to the problem of a single fictitious 
particle subjected to a perfectly reflecting barrier, as shown below. On the other hand, working on a lattice seems 
to be usually the most natural approach on physical grounds and is even a necessity when no continuous limit exists; 
an example of such a situation is the pure growth problem, equivalent to a directed walk, for which continuous limit 
of the master equation generates a purely mechanical Liouville equation, in which diffusive effects have disappeared 
Jl3| |. When the continuous limit exists, it is expected on physical grounds that both versions provide essentially the 
same results in the long time limit, when part of the microscopic details become irrelevant. Nevertheless, although 
leading terms are expected to coincide, subdominant corrections may play a rather important role, as shown below, 
since they usually follow power-laws in time with small exponents entailing that corrections are long-lived. When the 
relevant experimental timescale is short, results obtained by scaling hand-waving asymptotic arguments are of little 
physical interest and continuous models may even display serious shortcomings : as an example, as shown below, 
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the velocity at short times turns out to be infinite in the continuous approximation whereas its lattice analogue is 
perfectly well defined and is finite. In addition, somewhat surprinsingly, lattice effects persist even in the final regime, 
as constrasted to the ordinary diffusion of a single particle. This is why, as a whole, lattice models in continuous time 
are worthy to investigate: the discreteness of space, appearent in the transient as well as in the long-time dynamics, is 
not a minor feature of the problem. To be sure, lattice problems are by nature less "universal" than continuous ones 
in the sense that most results obtained in such a framework usually depend on microscopic details such as the lattice 
structure ( fll4|); nevertheless, they usually contain much more physically relevant information than continuous ones 
and, as such, can suggest high-resolved in time new experiments. 



II. THE MODEL 



The basic assumptions of the present lattice model are as follows : i) at some time (t = 0), a pair of particles is 
located on two given adjacent lattice sites labelled n = and 1. ii) when separated by more that one lattice spacing, 
each member of the pair has a symmetric diffusive motion, independently of the other; for simplicity, it is assumed 
that hopping can occur between one site and its nearest-neighbours. The hopping probability per unit time is denoted 
as W and allows to define a diffusion constant D = a 2 W, a being the lattice spacing and W~ x the diffusion time, iii) 
when the two particles arc located on two adjacent sites, each of them can only move onto the empty available site: 
this models the contact interaction. As a consequence, the two particles cannot stand on the same site and cannot 
cross each other. Any initial condition thus definitely induces a left-right symmetry breaking. Conventionally, with 
the above initial condition, the particle located at n — (resp. n = 1) at time t = will be given the label 1 (resp. 
2). 

With these assumptions, the master equation (M. E.) for the probability p n , m (t) to find the two particles 
located on sites n and m at time t can be written following the standard procedure. One first considers the evolution 
of p n ,m{t) between t and t + At, where At is a finite time interval, by exhausting all basic elementary jumps and 
expressing their probabilities in terms of the products W^At; by forming the difference p n . m (i + At) — p n ,m(t), by 
dividing both sides by At and taking the limit At — » 0, one eventually obtains the time-continuous following M. E. 
(Sn,m is the Kronecker symbol): 

Pn,m = - 2W [2 - (5 n ,m+l + S n ,m-l)] (1 - &n,m)Pri, m + W (p n +l,m +Pn,m-l)(l ~ ^n.m-l) (1 ~ ^n,m) 

+ W (p n -x,m + Pn,m+l){l - 5 n ,m+l) 0- - 8 n , m ) ■ (2.1) 

In such a way, as usual, all elementary steps involving two or more steps disappear (their probabilities scale as (At) fe 
with k > 2). The (1 — S nm ) factors account for the fact that the two particles cannot stand on the same site: when 
Pnn(O) = - which is the case fo r th e chosen initial condition -, then p nn (t) remains equal to zero at all times. 

Our aim is to solve eq. fl2.f| ) in order to deduce the full probability distribution {p n ,m(t)} and the related 
one-body (marginal, reduced) distributions {p n (l,t)} and {p m (2,t)} given by: 

+00 +00 

Pn(l,t) = E Pn,m(t) Pm(2,t) = ^ Pn,m(t) ■ (2-2) 



Eq. ( 2T ) can be solved by first introducing the generating (characteristic) function f(cj),ip,t) defined as: 



+30 



/(6V,*) - E Pn, m (t)^e^ . (2.3) 



n.m— — 00 



allowing to find each p n , m by inverse Fourier transformation or to get the marginal distributions by setting one of the 
arguments <f> or ifi equal to zero, and to find directly all the moments by successive derivations. It is readily seen that 
the characteristic function / satisfies the following homogeneous integro-differential equation: 

d f 2 * dd>' 

-/(0,V,O = -2VK(2-cos0-cosV>)/(0>,t)+2VF / [2 - cos0' - cos(0 + yj - (f>')}f(cj>' , cp + </> - <j>' , t) 



+ 2W — [2- cos(0- (/)')- cos (/.-cos (/) + ip - <j)' ,t) . (2.4) 







The above stated initial condition writes f(<p, ip, t = 0) = c 1 ^ and entails that p n ,m(t) — for n > m. By subsequently 
making a Laplace transformation: 
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r+oo 

f h (4>,^,z) = / dte- zt f(0,ip,t) , (2.5) 
Jo 

it is found that the function /l obeys the following homogeneous integral equation: 

/■2-7T J If 

zf L {cf>,i},z) - = -2W(2-cob0-cobV)M^^i*)+2W(2-cos0-oos^) / , (f> + ^ - 4>' , z) 

Jo 2vr 

+ W(sin^-sin^ cos^) /** ^ smtf htf ,</> + 1/> -</>', z) 

+ W(cos0-cos 2 ^A / 7^ cost' ,<j> + 4>-<j>', z) . (2.6) 
2 Jo 27T 

This is a Frcdholm equation with a degenerate kernel; as such, its solution is a linear combination of the three auxiliary 
quantities M(6), X{6) and Y{6) (6 = (/> + ip) defined as: 

M{9) = ^ ^/ L (0',0 + </;-0',z) 



2tt ji/ 



X{6) = Jo ^ W/l(<^ + ^-<M (2.7) 
= ^ 7^ sin07L(0'> + ^-0'^) 

M, X and Y" satisfy an inhomogeneous system which can be written down after a somewhat lengthy but straightforward 
calculation: 

(1 - cu)M - c 12 X - c 13 Y = cio 
-c 21 M + (I - c 22 )X ~ c 23 Y = c 2Q (2.8) 

- c 3 iM - c 32 x + (1 - c 33 )r = c 30 

where the various coefficients are equal to: 

e W/2 e -u 

C 10 — 7777 77i ;n\ — ~i — c 20 — 77777"; , n 777; i — r - ; — C30 



4VKcos(0/2) sinhu 8W | cos(0/2) | sinhu 8iW | cos(0/2) | sinhu 

1 - e-« I cos(0/2) I e"«- | cos(0/2) | tan(0/2) „ 

Cn = 1 taio\ 1 ■ £ 012 = ■ , C13 = .V e 11 - cos (0/2) 

I cos(0/2) I smhu smhu smhu 

e"" „ , t n 1 n\ 1 \ „ sinhu + cos 2 (6»/2) (e^ 11 - I cos(0/2) I) 

C21 = -v-r- 1 - coshu cos 0/2 c 22 = e"" 1 / / . 1 2.9 

sinliw I cos(0/2) I sinhu 

e- u sin(0/2) „ ...... 

C23 = c 32 = 77^-7- 1 \' l0 . 1 e cos (0/2) 

2 smhu I cos(0/2) | 

e-"tan(0/2) . _ n sinh u + sin 2 (0/2) (e""- | cos(0/2) |) 

C31 = . , 1 - coshu cos 0/2 c 33 = e \ . 1 — 

smhu I cos(0/2) | smhu 

In these equations, the quantity u is defined as follows: 

z + 4W , 

coshu = ——. ,„ 7-: 1 . (2.10 

4W I cos(0/2) I v ; 

u is a function of z and 0, uniquely defined by continuity from the branch assuming real positive values when z is a 
real positive number. 

The solution of the above system allows eventually to write down the Laplace transform of the generating 
function (2.5) as the following (Z = z/4W) : 



,. 1 e^cosgjH R{<t> + V>, Z) e^-V)/ 2 - | cos ±tH> , 

■ Ll '-'-'---- Z+l-i(cos0 + cosV) + V, Z) - cos 2 ±& """ 
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where: 



R(9, Z) = Z + 1 - [{Z + l) 2 - cos 2 (6>/2)] 



1/2 



(2.12) 



This is a first writing of the full solution of the present problem, since it allows to find all the moments of the two- 
body distribution {p n ,m(t)} in their Laplace representation. In addition, Laplace inversion of eq. (2.11) gives the 
characteristic function in the time domain: 



# /3(l-a)e^)/ 2 +/3 2 -a p4(a _ 1)m 



(3 2 - 2a + 1 
+ e # (1- p jW-^I^V 



-^e^j (^-x 2 ) 1 / 2 
3 7T x — a /3 2 — 2.t + 1 



In eq. (|2.13|) , T 3 denotes the Cauchy principal part and 



<b + ib (b — ib 
a = cos cos 



(3 = | cos 



(b + tp 



(2.13) 



(2.14) 



Moreover, by setting ip = (resp <j> = 0) in eq. ( [2.11 ) one readily obtains the Laplace transform of the characteristic 
function for the reduced (marginal) one-body distribution for the particle labelled 1 (resp. 2), fi^(<p,z) (resp. 
/2,l("0j z ))] their expressions for — 7r < <f> < +ir writes: 



1,L 



; | R(4>, Z) e^' 2 



Z + sin 2 | R(<b, Z) - cos 2 



(2.15) 



/ 2 ,l(^M) = ^[/u^z)] 



(2.16) 



The complex conjugation expresses the bias induced by the initial condition and simply means that the particle 1 
t ends , on the average, to move in one direction and the particle 2 in the opposite one. A first derivation at <p = of 
(2.15) gives the Laplace transform of the averaged coordinate of the particle 1, when the other can be found anywhere 
else: 



A second derivation gives the second moment: 



a 2 r %W 



814^ 



1/2' 



1/2 



8W 

z 



1/2 



- 1 



(2.17) 



(2.18) 



Inverse Laplace transformation now yields the exact expressions of the two first moments for the "left" particle (the 
one labelled 1) : 



< Xl > (t) = J {1 - c- T [(1 + 2T)7 (T) + 2Th(T)]} 



(2.19) 



< x 2 > (t) — a - < xi > (t) 



(2.20) 



<x 2 >(t) = %-T+ - {l-e- T [(1 + 2T)7 (T) + 2T7 1 (T)]} 



(2.21) 



In these equations, T is the dimensionless time T = 4Wt, whereas the I's denote the modified Bessel functions. The 
two mean square dispersions coincide and can be written in the symmetric form : 



Axl(t) = Axl(t) = 2a 2 Wt+ < x x > (t) < x 2 > (t) . 
It results that the mean distance between the two particles, d(t) = a — 2 < x\ > (t), is given by 



(2.22) 
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d(t) 



{l + c- T [(1 + 2T)I Q (T) + 2TI 1 (T)}} 



(2.23) 



At large times (Wt 3> 1), the preceding expressions have the following expansions : 

1/2 



< Xl > (t) 



-a I -Wt 

7T 



a 
4 



(2.24) 



Ax\{t) = Axj(t) = 2a 2 1 )Wt - a 



Wt 
27 



1/2 



3a? 
16 



(2.25) 



Equation ( 2.24 ) shows that, due to repeated collisions with the other particle, the average position of one particle 
eventually goes to infinity, although quite slowly (~ t 1 / 2 ) - and so does the average distan ce d(t): the repulsion 
induces an effective drift, characterized by a vanishing velocity at large times. Equation ( 2.25 ) displays the fact that 
the diffusion constant of each particle is reduced by the factor (1 — 1/t): the repulsive interaction partially inhibits 
the spreading, all the more since the first correction to the linear term is negative. 
On the opposite, for short times, t <C W -1 , one has : 



< Xl > (t) 



-aWt 



< x 2 > (t) ~ a (1 + Wt) 



(2.26) 



Ax\{t) = Ax 2 2 {t) 



a 2 Wt 



(2.27) 



Thus, at short times, the repulsion induces a normal drift with a finite velocity aW, whereas the diffusion 
constant is simply halved as compared to its value in the absence of interaction, since then each particle can essentially 
diffuse in half-space only, as in a directed walk. The two first moments in the lattice model are plotted in Fig. |l| 
(solid curves) for the "left" particle, using the exact expressions ( 2.1S| ), (2.21) and (2.22). 

The above results can be compared with those given by Sasamoto and Waditi || for the simplified version 
of the model in which each particle can move in one direction only (directed random walk). In such a case, each 
particle has obviously a finite velocity, as opposed to the present model where jumps can occur in both directions: 
the contact interaction induces for each particle a fluctuating boundary condition which gives rise to anomalous drift, 
the average coordinate increasing at all times as t 1 / 2 . More direct comparison can be done by considering the mean 
square dispersion. The quantities < X > and < X 2 > introduced by these authors (their eq. (4. 11)) are related to 
the above moments and are such that: 



< X 2 > - < X > 2 = -(Ax 
2 y 



Ax 2 ) + — 



(2.28) 



From (2.25), it is seen that, at large times and for the present model: 



< X 2 > - < X > 2 : 



2a \ I \Wt 

7T , 



(2.29) 



whereas the result of these authors for the directed walk reads: 



< X 2 > - < X > 2 : 



2a z 



(2.30) 



This shows that, as compared to the case without interaction, the diffusion coefficient is enhanced for the directed 
walk analyzed by Sasamoto and Waditi S whereas it is reduced for the both-way usual random walk model. 



III. THE CONTINUUM LIMIT 



The continuous space limit is worthy to analyze, although it cannot provide, by nature, the whole information 
contained in the lattice version and fails to represent transient regimes; the continuous space limit defines an oversim- 
plified framework and can only claim to describe features on large space and time scales. More precisely, if a is the 
order of magnitude of the underlying lattice and if r^iff denotes the diffusion time, wave vectors k and time-conjugate 
Laplace variables z are physically sensible only if they satisfy « a -1 and | z |<SC T diff- 
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FI G. 1. C ompa ri son o f the two first moments for the reduced one-bod y pro ba bility (left particle) in the lattice (solid lines, 
eqs. ( |2.19| ), ( 2.21 ), ( 2.22 )) and continuous framework (dashed lines, eqs. ( 3.1C ), ( 3.12 )). 



Generally speaking, the transition to continuous isotropic space is achieved from the square lattice framework 
by setting W — D /a 2 and by formally taking the limit a — > 0. Let us set = kia, ip = k 2 a; when the latter 
limit is taken, eq. (2.13) generates the Laplace transform of the characteristic function of the two-particle density, 
/ L (fci,/c 2 ,z): 



/ L (fci,fc 2 ,z) = 



1 



1 



[z + (D/2) (fci + k 2 y\ 1/2 [z + (D/2) (fci + k 2 f] 1/2 + i^Dj2 (k x - k 2 ) 



(3.1) 



This expression is much simpler than its analogous of the discrete version (comp. eq. ( 2.13 )). It can be reexpressed 
in terms of the total (center-of-mass momentum) K = k\ + k 2 and the relative (reduced) momentum k = {k\ — k 2 )/2 



fo{ki,k 2 ,z) 



1 



[z + (D/2)K 2 } 1/2 [z + (D/2)K 2 }' rz + iv^Dk 



21 1/2 



(3.2) 



displaying the fact that the diffusion constant for the center of mass motion is D/2, whereas the relative motion has 
a doubled diffusion constant 2D (indeed D plays the same role of an inverse mass in the reduction of a dynamical 
two-body problem). By using now the Efros theorem RDm, one easily obtains: 



f{h,k 2 ,t) = c 



= n -DK 2 t/2 



+°° d7 



e -r7(4t) e -iV2Dkr 



and Fourier inversion eventually yields the full two-body density probability P{x\,x 2l i) 



P(xi,x 2 ,t) = Q(x 2 -xi) 



1 



2nDt 



exp 



x^ -\~ x 2 
ADt 



(3.3) 



(3.4) 



where Q(x) is the Heaviside unit step function. This expression, obtained as the continuum limit of expression ( 2.15 ) 
is self-evident: provided that the inequality x\ < x 2 is satisfied, each particle has a free diffusion, independently of 
the other. Obviously, due to the (repulsive) interaction, P{x\,x 2 ,t) is not the product of two functions, one for each 
particle. Indeed, P(xi, x 2 ,t) is an integral of such products: 



P(xi,x 2 ,t) 



lim 



1 



dq 



eiq(X2 " xl) P -(*t + 4)mm 



e->o 2i7r 2 Dt J_ ao ' q - is 
exhibiting the obvious fact that the two particles are correlated at all times. 



(3.5) 



G 



The interparticle correlations are most simply expressed by the correlator C (t) , which is readily calculated from 
eq. (0: 



C(t) = < (x%- < xi >) (x 2 - < x 2 >) > = < xi > 2 



Dt 



(3.6) 



showing that the normalized correlator C(t)/ < Xi > 2 is constant in time: the correlations induced by the interaction 
never die out in one dimension. 
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FIG. 2. Reduced one-body density for the left particle in the continuous aproximation (see (3.7)). The diffusion constant D 
is taken as unity and each curve is labelled by the value of time; the asymmetry is less and less pronounced as time increases. 



be 



From cq. (3.4), the marginal probabilities in space-time for one particle, obviously asymmetric, are found to 



fi,2M) 



+ OC 



dx' P(x,x',t) 



-x 2 /(4Dt) 



l±$(x/V4Dt) 



(3.7) 



where the + (resp — ) sign refers to the right (resp. left) particle and where <!> denotes the probability integral 113]. The 



left particle marginal density is plotted in Fig. 
functions of the marginal densities: 



fi(q,t) 



Setting k 2 — (resp. ki = 0) in eq. ( p.3| ) yields the characteristic 

(3.8) 



f+OO J 



h(q,t) = [fi (g,t)]* ■ 

These functions give all the moments by successive derivations at q = 0; thus, at all times: 



<xi> (t) = - < x 2 > (t) 



2 -Dt 

7T 



< x\ > (t) =<x\> (t) = 2Dt 



Axl(t) = Axj(t) = 2(1 )Dt 

7T 



(3.9) 

(3.10) 
(3.11) 
(3.12) 
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These expre ssion s themselves reveal the simplifications carried out by the continuous approximation (comp. (2.19), 
( [2.21 ) and ( 2.22 )); they also display the expected fact that the exa ct re sults at all times of the continuous version 
coincide with the le ading asym ptoti c terms of the lattice version (see (2.24) and (2.25)); indeed, once the dimensionless 



time T — 4Wt in ( 2. IE ) and ( 2.21 ) is replaced by AD /a , the limit a — > automatically picks out the first leading 



term in the asymptotic expansion of the Bessel functions. Note however that lattice effects persist even at very large 
times: plots given in Fig. [j] allow the comparison between the lattice results (solid curves) and the continuous ones 
(dashed curves). For the mean square dispersion, the relative "error" decreases slowly as (Wt)^ 1 / 2 : it still amounts to 
10% for Wt — 100. The lattice effects indeed go to zero at infinite times, but in such a slow manner that subdominant 
terms cannot be neglected on reasonably large times. 

On the opposite, with no surprise, the behaviours at short tim es are markedly different (see ( 2.26| ) and FigJ|); 
for instance, the finite initial drift velocity in the lattice model (see (2.26), which could be observed with high-time 
resolution experiments, is infinite in the continuous approximation (eq. (3.10)). 

Note that the distribution function for the left particle (3.7) satisfies the conservation equation: 



dP x {x,t) 



dJ{x,t) 



with: 



J{x,t) 



dt 



-D 



dx 



dPi{x,t) 



dx 



2G 2 D (x,t) 



(3.13) 



(3.14) 



where Go denotes the normal density with a diffusion constant D. The additional term to the current is negative for 
all it must clearly be. 




FIG. 3. Comparison at short times of the two first moments for the reduced one-body probability (left particle) in the lattice 
(solid lines), and in the continuous framework (dashed lines). 



Obviously enough, the results of the continuum approximation are related to the simpler problem of a single 
fictitious brownian particle in the presence of a perfectly reflecting wall. For the two particles 1 and 2, the diffusion 
equation writes: 



d , . ( d 2 d 



—P(xt,x 2 ,t) = D ( + ) P(xi,X2,t) Vxi ^ x 2 , (3.15) 



the solution here satisfying the initial and boundary conditions P(xi,x 2 ,t — 0) = S(xi + 0)S(x2 — 0), 
P(%i,X2, t) =0 Vxi > X2- An easy way to solve such an equation is to transform to the center-of-mass frame 
by setting x — X2 — x±, X = \{x\ + x 2 ), P(xi, X2, t) — U(x,X,t); indeed, the diffusion constant D pla ys the 
same role as an inverse mass in a mechanical two-body problem. This allows to rewrite the diffusion equation ( 3.1 5| ) 
as the following : 
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d_ 

df 



U(x,X,t) = 



D d 2 
YdX 2 



2D 



dx 2 



U(x,X,t) ~ix 



(3.16) 



with U(x, X,t) = Vx < 0, VX £ R. 

Equation ( 0.16 ) displays the diffusion of the center-of-mass (with a constant D/2), decoupled from the diffusion 
of the relative coordinate (with a constant 2D) in the presence of a perfectly reflecting barrier at the origin. This 
means that the problem can be solved in two steps, by considering first the diffusive motion of a single fictitious 
particle with a doubled diffusion constant and by subsequently adding the effect of the center-of-mass motion. As a 
whole, the solution of the present problem in its continuous version is given by : 



P(xi,x 2 ,t) = G D/2 [ X = -{x\ +x 2 ),t ) P 2D (x =xi- x 2 ,t) 



(3.17) 



where Gp i 2 is the normalized gaussian function describing the free diffusive motion of the center of mass with the 
diffusion constant D/2, whereas P 2 d is the solution for the single particle (fictitious) constrained by the reflecting 
wall. As well known (Gardiner, 1990), if the fictitious particle stands at x = xo > with probability one at t = 0, 
one has: 



p 2D {x,t) = e(x) 



,+cc 2 
Jo * 



dk cos kxg cos kx e 



-IDhTt 



(3.18) 



By taking xq — 0+, by using (3.17) and by performing all the integrations, one recovers the expression given in (p.4| ), 
directly obtained as the limit taken from the lattice version. 

Note however that the present problem cannot be considered as a single particle one: the separation of the 
center-of-mass motion is trivial in the sense that it happens as the consequence of discrete space homogeneity, but 
this center-of-mass, itself having a diffusive motion (not a purely kinematical one), does participate to the sprea ding 
in such a way that correlations are always present; this is the reason why P(xi, x 2 ,t) does not factor out (see (3.4) 
and (3.5)) and sustains never-ending correlations, as simply expressed by (3.6). On a d ee per level, it is seen that 
the stochastic motion of each particle is no more Markoffian: the Fourier transforms ( |3.8| ), ( |3.9|) are not of the form 
q- a \i\ anc ; d G no t satisfy the Bachelier - Smoluchovski - Chapman - Kolmogorov chain equation. The presence of 
one particle thus strongly alters the nature of the motion of the other. 



IV. CONCLUSIONS 



The simple problem of two particles diffusing on a lattice with contact repulsive interaction was fully solved at 
all times by finding the two-particle characteristic function, from which marginal one-particle densities and moments 
can be easily derived. At short times, the repulsion induces a linear increase of the position for each particle, associated 
to a finite velocity, and a linear variation of the second cumulant, with a diffusion constant which is halved for obvious 
physical reasons. In the final regime, each particle still has an effective anomalous drift, its coordinate growing like 
i 1 / 2 , associated to a vanishing velocity. On the other hand, the mean square dispersion asymptotically follows a 
normal diffusion law; nevertheless, by comparison with its value in the absence of interaction, the diffusion constant 
is reduced by the non-trivial factor 1 — 1/n. As a whole, there is a cross-over between two normal diffusion regimes, 
the diffusion constant varying from D/2 at short times to (1 — 1/tt)D at large times. Independently of the long-lived 
subdominant terms, it is seen that each particle does not recover its "free" diffusion constant, even at infinite times. 
For sure, this is a low-dimensionnality effect. 

Taking the limit of continuous space first allows to make explicit connection with the problem of a fictitious 
particle in the presence of a totally reflecting barrier. Second, by analyzing this limit, it can be seen that marked 
differences exist between the lattice and the continuous versions. With no surprise, behaviours at short times are 
very different (the initial drift is finite in the lattice, infinite in the continuum). In addition, due to weakly decreasing 
subdominant terms, behaviours still substantially differ in both frameworks, even at large times: it can be said that 
lattice effects persist for physically relevant times, formally vanishing only at infinite times. Thus, when it is a given 
physical feature, the discreteness of space cannot be ignored. 
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Figure Captions 

Compariso n of t he two fi rst m oments for the reduced one-body probability (le ft pa rt icle) in the lattice (solid 
lines, eqs. ( 2.19D , (2.21), (2.22)) and continuous framework (dashed lines, eqs. ( |3.1C| ), ( 3 . 1 2| ) ) . 



2. Reduced one-body density for the left particle in the continuous aproximation (see (3.7)). The diffusion constant 
D is taken as unity and each curve is labelled by the value of time; the asymmetry is less and less pronounced 
as time increases. 

3. Comparison at short times of the two first moments for the reduced one-body probability (left particle) in the 
lattice (solid lines), and in the continuous framework (dashed lines). 
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